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Abstract 

We have studied the electronic and magnetic structures of the ternary iron arsenides AFe2As2 
(A = Ba, Ca, or Sr) using the first-principles density functional theory. The ground states of 
these compounds are in a collinear antiferromagnetic order, resulting from the interplay between 
the nearest and the next-nearest neighbor superexchange antiferromagnetic interactions bridged 
by As 4p orbitals. The correction from the spin-orbit interaction to the electronic band structure 
is given. The pressure can reduce dramatically the magnetic moment and diminish the collinear 
antiferromagnetic order. Based on the calculations, we propose that the low energy dynamics of 
these materials is described effectively by a t — Jh — J\ — ^2-type model 
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I. INTRODUCTION 



The recent discovery of superconductivity in LaFeAsO by partial substitution of O with F 
atoms below 26K[2] has stimulated great interest on the investigation of physical properties of 
iron-based pnictides. This type of quaternary compounds consists of alternative tetrahedral 
FeAs and LaO layers along the c-axis. The LaO layers act mainly as a charge reservoir. 
The superconducting pairing occurs in the FeAs layers. More recently, it was reported that 
the ternary iron-based arsenides AFe2As2 (A=Ba, Ca, or Sr) become superconducting upon 



hole or electron doping{3-6|. Similar as in LaFeAsO [jj, |8|, these ternary iron arsenides also 
exhibit a spin-density-wave-like anomaly and a structural transition from the tetragonal 
IA/mmm to the orthorhombic Fmmm group at some temperature between 140 K and 200 
K Furthermore, it was found that the high pressures can drive these undoped ternary 



iron arsenides superconducting 
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To investigate the mechanism of superconductivity in these materials, it is commonly 
believed that one needs to understand first the electronic and magnetic structures of the 
parent compounds. It has been shown that there is an essential similarity of electronic 
states nearby the Fermi level in AFe2As2 (A=Ba, Ca, or Sr) and LaFeAsO by studying 
the nonmagnetic state . Moreover, the density of states at the Fermi energy is only 
weakly doping dependent and the main effect of doping is a change in the relative sizes of 
the electron and hole Fermi surfaces 15j . 



In this paper, we report the electronic structures and magnetic orders and properties of 
AFe2As2 (A=Ba, Ca, or Sr) obtained from the first-principles electronic structure calcula- 
tions. By comparison of the energy of the non-magnetic state with those of a number of 
magnetic ordered states, we find that, similar as in LaFeAsO, the ground state of AFe2As2 
is in a collinear antiferromagnetic order. We have also studied the spin-orbit interaction and 
the pressure effect in these materials. The electronic and magnetic structure is found to be 
strongly affected by the pressure effect, but weakly by the spin-orbit interaction. 

II. COMPUTATIONAL APPROACH 

AFe2As2 (A=Ba, Ca, or Sr) takes the ThCr2Si2 type structure and an "A" layer plays a 
similar role as a LaO layer in LaFeAsO. Although AFe2As2 can be considered as a tetragonal 
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FIG. 1: (Color online) AFe2As2 (A=Ba, Ca, or Sr): a) a tetragonal unit cell containing two formula 
units; b) a primitive unit cell containing one formula unit; c) schematic top view of the FeAs layer. 
The small red dashed square is an a x a unit cell, while the large black dashed square is a \/2a x y/2a 
unit cell. The Fe spins in the collinear antiferromagnetic order are shown by arrows. 

crystal with two formula units included in the corresponding unit cell as shown in FigJTfa), 
its primitive unit cell is constructed by considering AFe2As2 as a triclinic crystal, in which 
only one formula unit cell is included as shown in Fig. Hfb). In the calculations, we adopted 
the primitive cell as calculation cell with the experimental lattice constants as the input 
parameters. In the calculation of the electronic structures of the nonmagnetic, the ferro- 
magnetic and the square antiferromagnetic Neel states, the a x a FeAs cell is taken as the 
base cell (shown in Fig. QJc)). In the calculation of the collinear antiferromagnetic state, 
the unit cell is doubled and the base cell is the v2a x V2a FeAs cell as shown in Fig. [TJc). 



In our calculations the plane wave basis method was used 



gradient approximation (GGA) of Perdew-Burke-Ernzerhof 17| for the exchange-correlation 



161 ] . We used the generalized 



potentials. The ultrasoft pseudopotentials 



18] were used to model the electron-ion inter- 



actions. After the full convergence test, the kinetic energy cut-off and the charge density 
cut-off of the plane wave basis were chosen to be 600eV and 4800eV, respectively. The Gaus- 
sian broadening technique was used and a mesh of 16 x 16 x 8 k-points were sampled for the 
Brillouin-zone integration. The internal atomic coordinates within a cell were determined 
by the energy minimization. 



3 



III. RESULTS AND ANALYSIS 



A. Nonmagnetic State 

We first studied the nonmagnetic state of the compound AFe 2 As 2 (A = Ba, Sr, or Ca), 
which is the high temperature phase of these materials. The electronic band structure of 
this state also provides a reference for studying of the low temperature magnetic phases. 
This can help us to understand the mechanism or the interactions that drive the magnetic 
phase transition and the related structural transition [lqI ]. 



1. BaFe 2 As2 

In our calculation, the experimental tetragonal crystal lattice parameters a = b = 3.9625A 



and c = 13.0168A [9] were adopted for BaFe 2 As 2 . Figs. |2] and [3] show the calculated 
density of states (DOS) of BaFe 2 As 2 in the nonmagnetic, square antiferromagnetic Neel, 
and collinear antiferromagnetic states, respectively. As revealed by Fig. |2(a), similar to 
LaFeAsO, the density of states of BaFe 2 As 2 consists of mainly the Fe-3<i states from -2eV to 
2eV around the Fermi energy. Further analysis of the calculation (Fig. [3^) shows that the 
crystal field splitting of the Fe-3d orbitals is much smaller than the one in transition metal 
oxides. This is not very surprising since the electronegativity of As is much weaker than O. 
This suggests that the Fermi surface may have the contribution from all the Fe-3d orbitals. 

The electronic band structure and the Fermi surface of BaFe 2 As 2 in the nonmagnetic 
state are shown in Fig. HI As shown in Fig. H^a), there are four Fermi surface sheets, 
contributed from the four bands crossing the Fermi energy (Fig. H](a)). Among them, the 
two cylinder-like Fermi-surface sheets centered around X-P are from the electron bands. 
They correspond to the two Fermi surface sheets around M-A in LaFeAsO. The other two 
Fermi surface sheets centered around T-Z are from the hole bands. These results agree 



qualitatively with the experimental observation 



20] 



The e nergy d ispersion of the electronic bands along the c-axis is much larger than that in 



LaFeAsO 
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22( | . The sectional views parallel to (001) plane through Z point and T point are 



different due to the large dispersion (Figs. Etc) and (d)). There is one band just below the 
Fermi energy along T-Z in Fig. H(a), which corresponds to the third hole Fermi surface sheet 
given m Ref . Q. The volumes enclosed by these Fermi surface sheets are 0.26 electrons/cell 




FIG. 2: (Color online) Total and orbital-resolved partial density of states of BaFe2As2 per formula 
unit in a) nonmagnetic state; b) square antiferromagnetic Neel state with spin-up (the same for 
the spin-down); and c) collinear antiferromagnetic state with spin- up (the same for the spin-down), 
respectively. 

and 0.26 holes/cell, respectively. The electron carrier concentration is the same as the hole 
carrier concentration. Both are equal to 2.54 x 10 21 / cm 3 . The compound BaFe2As2 is thus 
a semimetal with a low carrier concentration between normal metals and semiconductors, 



similar to what we found in LaFeAsO 22]. The density of states at the Fermi energy is 3.93 
state per eV per formula unit cell. The corresponding electronic specific heat coefficient 7 
= 9.26mJ/(K 2 * mol) and the Pauli paramagnetic susceptibility \v = 1-60 x 10 _9 m 3 /mol. 
These calculated physical quantities are also summarized in Table HJ which are well close to 
the experimental values|9, 

We study the plasma excitation in the semimetal BaFe2As2 as well. The plasma frequency 
f p is computed as follows, 



27r y m h m* 

where nn and n e are the hole and electron carrier densities, respectively and m* h and m* are 
the effective masses of hole and electron, respectively. We use the second order polynomial 
fitting to estimate the effective masses of carriers from the calculated band structures. We 
obtain the hole and electron effective masses as m* = 0.60m e and m* h = 0.86m e with m e 




FIG. 3: (Color online) Total and projected density of states at the five Fe-3cf orbitals of BaFe2As2 
per Fe atom in a) nonmagnetic state, b) square antiferromagnetic Neel state, and c) collinear 
antiferromagnetic state, respectively. 




FIG. 4: (Color online) (a) Electronic band structure of BaFe2As2 in a primitive unit cell for the 
nonmagnetic state. The Fermi surface (b) and its sectional views through symmetrical k-point Z 
(c) and T (d) parallel to (001) plane. 



being electron mass. Because of irregular and strongly anisotropic band structure, the 
estimated effective masses are with not small uncertainty. We then calculate the plasma 
frequency as f p = 25392 cm -1 in the nonmagnetic state. These values are also listed in 
Table [I] in comparison with the available experimental values. 
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TABLE I: Calculated (Cal.) physical quantities in comparison with the available experimental 
(Exp.) values for AFe2As2 (A=Ba, Sr, or Ca) in the nonmagnetic state (NM) and the collinear 
antiferromagnetic state (Col). The units for the carrier density p, the plasma frequency f p , the 
electronic specific heat coefficient 7, and the Pauli paramagnetic susceptibility Xp are 10 21 /cm 3 , 
cm" 1 , mJ/(K 2 *mol), and 10 _9 m 3 /mol, respectively. J\ and J2 (meV/S* 2 per Fe) are the superex- 
change antiferromagnetic couplings between the nearest and the next nearest neighbor Fe spins S in 
a6-plane, respectively while J z (meV/5 2 per Fe) is the interlayer superexchange antiferromagnetic 
coupling of Fe spins S perpendicular to a6-plane. 



AFe 2 As 2 


State 


P 


fp 
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Coupling 


hole 


electron 


Cal. 


Exp. 


Cal. 


Exp. 


h 


J2 


Jz 


BaFe2As2 


NM 


2.54 


2.54 


25392 


12900[27] 


9.26 


37[28j, 16[9j, 


1.60 


25.5 


33.8 


3.1 


Col 


0.10 


0.21 


7717 


4660 [27] 


5.68 




SrFe2As2 


NM 


3.33 


3.33 


27386 


13840[27] 


7.71 


6.1 [23] 
6.5 


1.33 


14.7 


33.4 


7.8 


Col 


0.13 


0.04 


4249 


4750 [27] 


3.63 




CaFe 2 As 2 


NM 


4.21 


4.21 


20883 




9.31 


8.2 


1.60 


-1.4 


26.8 


14.8 


Col 


0.09 


0.09 


5557 




2.62 





2. SrFe 2 As2 



Like BaFe2As2, SrFe2As2 can also become superconducting by chemical doping or by ap- 
plying a high pressure. These two compounds have the same crystal structure. The tetrago- 
nal crystal lattice constants determined by experimental measurements are a = b = 3. 9259 A 
and c = 12. 375 A for SrFe 2 As2 [4]. They are smaller than the corresponding parameters for 
BaFe2As2, especially along the c-axis. The electronic band structure and the Fermi surface 
shown in Fig. [5J are similar as the ones for BaFe2As2- The difference is that the band along 
T-Z, which is just below the Fermi energy in BaFe2As2, moves slightly upward and intersects 
with the Fermi level. This results in the third hole-type Fermi sheet centered around T-Z. 
Therefore, in SrFe2As2, there are three hole-type and two electron-type Fermi surface sheets. 

The volumes enclosed by these Fermi sheets give 0.32 electrons/cell and 0.32 holes/cell for 
SrFe2As2- The electron (or hole) carrier concentration is about 3.33 x 10 21 /cm 3 . The density 
of states at the Fermi energy is about 3.27 state per eV per formula unit. The corresponding 
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FIG. 5: (Color online) Calculated electronic structure of SrFe2As2 in a primitive unit cell for the 
nonmagnetic state, (a) Electronic band structure; The Fermi surface (b) and its sectional views 
through symmetrical k-point Z (c) and T (d) parallel to (001) plane. The symmetrical k-points in 
the Brillouin zone are referred to Fig. S^b). 

electronic specific heat coefficient and Pauli susceptibility are 7 = 7. 71m J/ (K 2 * mol) and 
X P = 1.33 x 1CT 9 'm 3 /mol, respectively. And the plasma frequency f p is computed as about 
27386 cmT 1 . These calculated quantities are also reported in Tabled! The total and Fe-3<i 
projected density of states of SrFe2As2 are similar as the ones for BaFe2As2. The low energy 
excitations are also dominated by Fe-3d orbitals from -2eV to 2eV around the Fermi energy. 



3. CaFe 2 As 2 



The tetragonal crystal lattice parameters obtained by experiments [261] . a = b = 3.912A 
and c = 11.667A, are used in our calculation for CaFe2As2- Fig. [6]shows the electronic band 
structure and the Fermi surface. Similar to SrFe2As2, there are three hole- type and two 
electron- type Fermi surface sheets for CaFe2As2. Because the third hole-type Fermi surface of 
CaFe 2 As2 expands into a cylinder-like shapes centered around T-Z, the cross section through 
T and X in (001) plane has one more cutting line than that of SrFe2As2. From the volumes 
enclosed by these Fermi sheets, we determine the electron and hole concentrations are 0.38 
electrons/cell and 0.38 holes/cell, respectively. The corresponding electron (or hole) carrier 
density is about 4.21 x 10 21 /cm 3 . The density of states at the Fermi energy is 3.95 state 
per eV per formula unit, and the electronic specific heat coefficient 7 = 9.31m J/ (K 2 * mol) 
and Pauli paramagnetic susceptibility Xp — 1-60 x 10 _9 m 3 /mol. And the plasma frequency 



8 




FIG. 6: (Color online) Calculated band structure and Fermi surface of CaFe2As2 in a primitive unit 
cell for the nonmagnetic state, (a) Electronic band structure; (b) Fermi surface; And the sectional 
views through symmetrical k-point Z (c) and T (d) parallel to (001) plane. The symmetrical 
k-points in the Brillouin zone are referred to Fig. EKb). 




FIG. 7: (Color online) The electronic structure of BaFe2As2 in a primitive unit cell for the 
metastable square antiferromagnetic Neel state, (a) Electronic band structure; The Fermi sur- 
face (b) and its the sectional views through symmetrical k-point Z (c) and T (d) parallel to (001) 
plane. 

f p is computed as about 20883 cm -1 . These calculated quantities are summarized in Table 
[H The total and Fe-3<i orbital projected density of states of CaFe2As2 are almost the same 
as the ones for SrFe2As2. 
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B. Antiferromagnetic Neel State 



To study the electronic structures of AFe2As2 in a magnetic state, we break the spin up- 
down symmetry by assigning a finite magnetic moment to each Fe atom. The ferromagnetic 
state is found to be not stable in AFe2As2, similar as in LaFeAsO 22J. However, the 



antiferromagnetic Neel state is meta-stable. Its energy is lower than the corresponding 
nonmagnetic state by 0.232 eV, 0.175 eV, and 0.133 eV per formula unit for BaFe2As2, 
SrFe2As2, and CaFe2As2, respectively. The magnetic moments in BaFe2As2, SrFe2As2, and 



CaFe 2 As 2 are found 
the one for LaFeAsO 



;o be respectively 2.3 /ig/Fe, 2.2 fj, B /Fe, and 2.0 fis/Fe, similar as 



221 ] . while a moment of about 1.7 /i^/Fe was found for all these three 



systems in the full potential local-density approximation calculations |29|. The density of 
states for the antiferromagnetic Neel state of BaFe2As2 is shown in Figs. E(b) and [3(b). 
The corresponding electronic band structure with the Fermi surface shapes is shown in Fig. 
[JJ The electronic structures and the Fermi surface topology are similar for the other two 
ternary iron arsenides in the antiferromagnetic Neel state. 

As shown in Fig. [5(b), the density of states around there is substantially reduced around 
-0.5eV in comparison with that in the nonmagnetic state (Fig. 12(a)). The missing states are 
pushed down to around -2.0eV. This change can be also seen by comparing the electronic 
band structures of the nonmagnetic and the square antiferromagnetic Neel states, shown in 
Fig. E(a) and Fig. [7(a). In contrast to the nonmagnetic state, there are only two bands 
crossing the Fermi energy, forming a ring-like hole-type shape around Z, a cylinder-like 
electron-type sheet centered around X-P, and a small-pocket-like electron-type sheet around 
Z. 



C. Collinear Antiferromagnetic State 

For BaFe2As2, SrFe2As2, and CaFe2As2, the antiferromagnetic Neel state is a meta-stable 
state among the rich magnetic structures. It is found that the true ground state is in fact 
a collinear antiferromagnetic state with the interlayer Fe moments in antiferromagnetic 
alignment. The spin configuration of this state in FeAs layer is schematically shown in Fig. 

Die). 
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1. BaFe 2 As2 



For BaFe2As2, the energy of the collinear antiferromagnetic state is lowered by 0.400 eV 
per formula unit than that of the nonmagnetic state. The magnetic moment is found to be 
about 2.65/ie for each Fe ion in this state. The electronic band structure and the density 
of states are shown in Fig. [B^a) and Fig. E(c), respectively. In contrast to the nonmagnetic 
state, the density of states of As Ap orbitals in the low energy range from -2eV to OeV is 
substantially enhanced in the collinear antiferromagnetic state. 

Fig. M^c) shows the Fermi surface shapes of BaFe2As2 in the collinear antiferromagnetic 
state. There are two bands crossing the Fermi level. They lead to a hole-type Fermi 
surface and an electron-type Fermi surface, respectively. Here we remind that if we consider 
BaFe2As2 by a tetragonal unit cell, the axis MA in the Brillouin zone of the tetragonal unit 
cell is now folded into the TZ, in which M point coincides with T point. This means that the 
electron-type Fermi sheet around M point in the tetragonal unit cell in the nonmagnetic state 
will be gapped when the collinear antiferromagnetic order takes place. From the volumes 
enclosed by these two Fermi sheets, we find that the electron carrier concentration is 0.042 
electrons/cell, namely 2.05 x 10 20 /cm 3 . The corresponding hole carrier concentration is 0.021 
holes/cell, namely 1.03 x 10 20 /cm 3 . Thus the electron concentration dominates over the hole 
concentration, unlike in the nonmagnetic state. By using the second order polynomial to fit 
the calculated band structures, we estimate the electron effective masses as m* = 0.31m e 
with m e being electron mass. We then calculate the plasma frequency as f p = 7717 cm -1 . 
Here the contribution of the holes to the plasma frequency is negligible. These values are 
also listed in Table |T] in comparison with the available experimental values. 

In comparison with the nonmagnetic state, the electron effective mass becomes much 
lighter. Meanwhile, the carrier density is reduced by an order of magnitude, as shown in 
Table [B This reduction is significantly smaller than that in LaFeAsO. In LaFeAsO, the 
carrier density is reduced by two orders of magnitude in the collinear antiferromagnetic 
state than in the non-magnetic state 19]. The density of states at the Fermi energy is 
2.41 state per eV per formula unit cell. From this, we find that the electronic specific heat 
coefficient 7 = 5.68mJ/(K 2 * mol). It should be emphasized that the calculated value of 
7 in the collinear antiferromagnetic state should be smaller than the intrinsic value of 7 at 
zero temperature measured by experiments since the low energy quantum spin fluctuations 
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FIG. 8: (Color online) Calculated band structure and Fermi surface of BaFe2As2 in the collinear 
antiferromagnetic state with the antiparallel alignment between the interlayer Fe moments along 



c-axis 
surface. 



30l | . which is the ground state, (a) Electronic band structure; (b) Brillouin zone; (c) Fermi 



are suppressed in the calculations, as discussed in Section |TV] Therefore, the measured 7 in 
the low temperature limit should be bound between the two calculated 7 values obtained 
by our calculations for the nonmagnetic and collinear antiferromagntic states, respectively. 
In Table HI we compare the calculated specific heat coefficients 7 with the corresponding 
experimental values. For BaFe2As2, there is a large variety in the measurement values. 
Experimentally, the first reported measurement value of 7 is about 37 mJ/(K 2 * mol) 28], 
much larger than our DFT result. However, with the sample quality being improved, the 
second reported value of 7 becomes 16mJ/(K 2 * mol) and the latest reported value of 
7 is 6.1 mJ/(K 2 * mol) [23], which is in fact good consistent with our DFT calculation. For 
SrFe 2 As 2 and CaFe 2 As 2 , our calculated results are also consistent with the measurement 
values. 

In the collinear antiferromagnetic state, there is further a small energy gain if the Fe- 
Fe distance is reduced along the spin parallel alignment direction and expanded along the 
spin anti-parallel alignment direction. This leads to a structural transition from tetragonal 
space group IA/mmm to orthorhombic space group Fmmm, similar to that in LaFeAsO 



191 ] and observed by the neutron scattering 10]. This lattice relaxation is energetically 



favorable because the direct ferromagnetic exchange favors a shorter Fe-Fe separation while 
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FIG. 9: (Color online) Calculated band structure and Fermi surface of SrFe2As2 in the collinear 
antiferromagnetic state with the antiparallel alignment between the interlayer Fe moments along 



c-axis 



30|. (a) Electronic band structure; (b) Side view: Fermi surface; (c)Top view: Fermi surface. 



The symmetrical k-points in the Brillouin zone are referred to Fig. E^b). 

the antiferromagnetic superexchange favors a larger Fe-As-Fe angle. It turns out that the 
angle 7 in ab-plane is no longer rectangular (Fig. [He)) and the energy gain is 3 meV per 
formula unit, when 7 is about 90.8°, similar as in LaFeAsO [19]. The correction from this 
lattice distortion to the electronic band structure as well as the Fe moments is very small. 

Along the c-axis, we find that the Fe-spins between the nearest neighbor FeAs layers 
interact antiferromagnetically and the energy gain by taking the anti-parallel alignment 
is about 0.012eV per formula unit cell in comparison with the parallel alignment. This 
antiferromagnetic interaction between the nearest neighboring FeAs layers is significantly 
larger than that in LaFeAsO 19]. In Ref. Q, we reported the electronic band structure and 
the Fermi surface of AFe2As2 (A=Ba, Sr, Ca) in the collinear antiferromagnetic order with 
the parallel alignment along c-axis. 



2. SrFe 2 As 2 

For SrFe2As2, the energy of the collinear antiferromagnetic state is lowered by 0.383 eV 
per formula unit than the one of the nonmagnetic state. The magnetic moment is about 
2.55/zb per Fe atom. The electronic band structure and the Fermi surface are shown in Fig. 
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FIG. 10: (Color online) Calculated band structure and Fermi surface of CaFe2As2 in the collinear 
antiferromagnetic state with the antiparallel alignment between the interlayer Fe moments along 



c-axis 



301 ] . (a) Electronic band structure; (b) Side view: Fermi suface; (c)Top view: Fermi surface. 



The symmetrical k-points in the Brillouin zone are referred to Fig. 0(b) . 



[9j From the volumes enclosed by the Fermi surface sheets, we find that the electron and 
hole carrier concentrations are about 0.008 electrons/cell and 0.025 holes/cell, namely about 
0.44 x 10 20 /cm 3 and 1.31 x 10 20 /cm 3 , respectively. It turns out that the hole concentration 
dominates over the electron concentration, unlike the case of BaFe2As2. And the plasma 
frequency f p is computed as about 4249 cm -1 . In comparison with the nonmagnetic state, 
the carrier density is much reduced, by more than an order of magnitude. The density of 
states at the Fermi energy is 1.54 state per eV per formula unit cell, and the electronic 
specific heat coefficient 7 = 3. 63m J / (K 2 *mol). These calculated quantities are also reported 
in Table HI 



3. CaFe 2 As 2 

For CaFe2As2, the collinear antiferromagnetic ordering can lower the ground state energy 
by 0.352 eV per formula unit in comparison with the nonmagnetic state. There are about 
2.49/is moment around each Fe atom. The electronic band structure and the Fermi surface 
are shown in Fig. [TUJ The electron (hole) carrier density is about 0.016 electrons/cell (0.016 
holes/cell), or 0.90 x 10 20 /cm 3 (0.90 x 10 20 /cm 3 ). Thus the electron concentration and the 
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FIG. 11: (Color online) Electronic band structures of BaFe2As2 including spin-orbit interaction in 
the nonmagnetic state. 

hole concentration are in balance, different from the ones in BaFe2As2 and SrFe2As2. And 
the plasma frequency f p is computed as about 5557 cm -1 . Similar to SrFe 2 As 2 , the carrier 
density is also much reduced by more than an order of magnitude in comparison with the 
nonmagnetic state. The density of states at the Fermi energy is about 1.11 state per 
eV per formula unit cell and the electronic specific heat coefficient 7 = 2.62mJ/(K 2 * mot). 
These calculated quantities are also reported in Table [fl 

D. Spin-Orbit Interaction 

The spin-orbit interaction results from the relativistic effect. It is known that this inter- 
action leads to many interesting features in transition metal oxides. To study how strong 
this interaction can affect electronic properties of these materials, we performed a relativis- 
tic calculation for BaFe2As2- Fig. [TT1 shows the electronic band structures of BaFe2As2 by 
including spin-orbit interaction in the nonmagnetic state. By comparison with Figs. H]^a) 
and[S](a), we find that the spin-orbit interaction splits the band mainly along T-Z into two 
bands by about 50~150 meV around the Fermi energy. 
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FIG. 12: (Color online) a) Calculated energy differences of BaFe2As2 between the magnetic state 
and the nonmagnetic state (NM) per formula unit cell under various pressures. Note the energy of 
the nonmagnetic state without pressure is set to zero, b) The magnetic moments for the collinear 
antiferromagnetic state (Collinear) and the antiferromagnetic state (AFM) under the pressures. 



E. Pressure Effect 



Experimentally it was reported 
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lMl3l| that the ternary iron arsenides AFe2As2 (A=Ba, 
Ca, or Sr) can become superconducting under high pressures without doping. This provides 
another route to study the salient features of these materials in connection with the intrinsic 
electronic structures and properties. The superconducting phase appears when the pressure 
is in the ranges of 2.5~8, 28~37, and 22~58 kbar for CaFe2As2, SrFe2As2, and BaFe2As2 re- 
spectively. The highest superconducting transition temperatures for these three compounds 
are 12K, 27K, and 29K, respectively. 

Fig- IT2T a) shows the pressure dependence of relative energies for the nonmagnetic, square 
antiferromagnetic Neel and collinear antiferromagnetic states of BaFe2As2, respectively. The 
collinear antiferromagnetic state is robust against the pressure and has the lowest energy 
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Note that 



in the pressure range studied, consistent with the result reported in in Ref. 
in this calculation, the lattice parameters are relaxed when the pressure is zero. Thus the 
energy differences are slightly different from that given before. 

The lattice structure of FeAs layers, including the angle between Fe and As atoms and 
the bond lengths of Fe-Fe, Fe-As, and As-As, is hardly changed by the pressure from to 
80 kbar. However, the contraction along the c axis is much more pronounced, especially in 
the nonmagnetic state. This indicates that the distance between the neighboring FeAs and 
Ba layers is substantially reduced by pressure. This can strengthen the coupling between 
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the neighboring FeAs layers and stabilize the long range antiferromagnetic correlation in the 
FeAs layer. 

The magnetic moment in the magnetic states decreases almost linearly with pressure, as 
shown in Fig. [L2T b). In our calculations, the magnetic moment is always larger than in 
the collinear antiferromagntic state. The superexchange interaction thus remains dominant 
when a pressure is imposed. 



IV. EFFECTIVE MODEL 



A. Local Moment versus Itinerant Electrons 



Physically the moments of Fe ions result from the on-site Coulomb repulsion and the 
Hund's rule coupling of 3d orbitals. An isolated Fe ion in a 2+ valency has a spin S — 2 
with a large magnetic moment of 4/xb. In iron-pnictide semimetals, the effective moment 
of a Fe ion (in a 2+ valency) will be reduced by its hybridization with other atoms and by 
the Coulomb screening of itinerant electrons. However, it will remain finite if the Hund's 
rule coupling and the on-site Coulomb repulsion is strong enough in comparison with the 
hybridization and other screening effects, as we found in LaFeAsO [h]]. 

In our calculations, by projecting the density of states onto the five 3d orbitals of Fe in the 
collinear antiferromagnetic state of BaFe2As2 (Fig. E(c)), we find that the five 3d orbitals 
of Fe are almost completely filled by up-spin electrons and nearly half-filled by down-spin 
electrons in one of the two sublattices (or completely filled by down-spin electrons and half- 
filled by up-spin electrons in the other sublattice). This indicates that the crystal field 
splitting imposed by As atoms is very small and the Fe 3c/-orbitals hybridize strongly with 
each other. We can see that this is a universal feature for all iron pnictides, as we first found 
in LaFeAsO [19(]. The strong polarization of Fe magnetic moments is thus due to the Hund's 
rule coupling. 

In low temperatures, the Fe moments will interact with each other to form an antifer- 
romagnetic ordered state. These ordered magnetic moments have been observed by elastic 
neutron scattering and other experiments [7, 1(3] . However, they are not exactly the mo- 
ments obtained by the DFT calculations, as we indicated first for LaFeAsO 19|. This is 
because the DFT calculation is done based on a small magnetic unit cell and the low-energy 
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quantum spin fluctuations as well as their interactions with itinerant electrons are frozen by 
the finite excitation gap due to the finite-size effect. Thus the moment obtained by the DFT 
is the bare moment of each Fe ion. It should be larger than the ordering moment measured 
by neutron scattering and other experiments. Our calculations show that the bare magnetic 
moment around each Fe atom is about 2.2 ~ 2.6 in all iron pnictides and in different 
magnetically ordered states. 

In high temperatures, there is no net static moment in the paramagnetic phase due to 
the thermal fluctuation, but the bare moment of each Fe ion can still be measured by a fast 
local probe like ESR (electron spin resonance). Very recently, the bare moment of Fe has 
been observed in the paramagnetic phase by the ESR measurement j^. The value of the 
moment detected by ESR is about 2.2 ~ 2.8 /xg i n good agreement with our DFT result, 
which is but significantly larger than the ordering moment in the antiferromagnetic phase. 



19j . from the spatial distributions of elec- 



Again similar to what we found in LaFeAsO 
trons, we further find that there is a strong hybridization between neighboring Fe and As 
ions. This strong hybridization can mediate an antiferromagnetic superexchange interac- 
tion between the Fe moments. This superexchange interaction is antiferromangtic since the 
intermediated state associated with the hopping bridged by As ions is a spin singlet. On 
the other hand, there is a relatively small but finite hybridization between two neighboring 
Fe ions. This direct hybridization of Fe 3d orbitals can induce a direct exchange interac- 
tion between the Fe moments. This direct exchange interaction is ferromagnetic due to 
the Hund's rule coupling. Furthermore, there is a strong covalent bonding or hybridization 
between As 4p orbitals although the separation between As atoms is relatively large. This 
hybridization gives rise to a broad As 4p band below the Fermi level. Thus As 4p states 
are not truly localized. They form an electron network connecting As ions through covalent 
bonding, similar to what we found in iron chalcogenides (a-FeSe and a-FeTe) 32j. However, 
there is difference that the band formed by As 4p-orbitals is insulating, as shown in Fig. 
[21 while the band formed by Te 5p-orbitals is metallic (see Ref. S). It turns out that the 
exchange interaction bridged by As 4p orbitals is short ranged (just J± and J 2) while the 
one bridged by Te 5p orbitals can be long ranged, as we found in a-FeTe in which there is 
a substantial third nearest neighbor antiferromagnetic superexchange interaction J 3 besides 
J\ and J2 [321 ] . 
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B. t — Jh — J\ — Ji Hamiltonian 



For iron-based pnictides, the low energy spin dynamics could be approximately described 
by an antiferromagnetic Heisenberg model with the nearest and the next-nearest neighbor 
exchange interactions. However, the Fe spin (or magnetic moment) is not quantized since 
the electrons constituting the moment can propagate on the lattice, like in hole or electron 
doped high T c superconductivity cuprates. Besides these superexchange interactions, the 
on-site Hund's rule coupling among different Fe 3c? orbitals is important. This is because 
the crystal splitting of Fe 3d levels is very small and the spins of Fe 3d electrons are polarized 
mainly by this interaction. Thus we believe that the low-energy physical properties of these 
iron-based pnictides can be approximately described by the following effective Hamiltonian 

H — tffclaCjp - J H ^2 ^ia ■ Si/3 

{ij),a/3 hoe£fl 

+ Jl Si a ■ Sjp + J2 Si a ■ Sjp, (2) 

where (ij) and ((ij)) represent the summation over the nearest and the next-nearest neigh- 
bors, respectively, a and are the indices of Fe 3d orbitals. c\ a is the electron creation 
(annihilation) operator. 

Sia — c ia^ C io! (3) 

is the spin operator of the a orbital at site i. The total spin operator at site i is defined by 
Si = J2 a Sia- In Eq. (j2J), Jh is the on-site Hund's coupling among the five Fe 3d orbitals. 
The value of Jh is generally believed to be about 1 eV. are the effective hopping integrals 



33|. 



that can be determined from the electronic band structure in the nonmagnetic state 

The nearest and the next-nearest neighbor antiferromagnetic coupling constants, J\ and 
J2, in Eq. (j2j) can be calculated from the relative energies of the ferromagnetic, square 
antiferromagnetic, and collinear antiferromagnetic states with respect to the non-magnetic 
state. For the corresponding detailed calculations, please refer to the appendix in our paper 
111 Ref. [l9. For BaFe2As2, the energy of the ferromagnetic state is about 8meV per Fe 
lower than the nonmagnetic state and the Fe moment is about 2.3/z#. From this and the 
relative energies of the two antiferromagnetic states, we find that the exchange constants 
are approximately given by J\— 25.5 meV / S 2 per Fe and J 2 = 33.8 meV/S 2 per Fe (S is the 
spin of the Fe ion) for BaFe2As2. Meanwhile, we also find that the interlayer superexchange 
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antiferromagnetic coupling J z =3.1 meV/S 2 per Fe. In obtaining these values, we have 
assumed that the contribution of itinerant electrons to the energy is almost unchanged 
in different magnetically ordered states. Since the bare and Jh can be considering 
independent of magnetic structures, the relative energies between different magnetic states 
are not affected by itinerant electrons. For SrFe2As2 and CaFe2As2, the values of J\ and 
J2 with J z are determined as well and given in Table [B As we notice, J\ in CaFe2As2 is 
negative, namely ferromagnetic, which is unique among iron pnictides. 

C. Discussion 

We plot the superconducting critical temperatures T c versus J2 for these compounds with 
doping carriers or by applying high pressures in Fig. [13J Interesting, as we see, the maxi- 
mum critical temperatures T c are in proportion to the next-nearest neighbor superexchange 
interaction J2. This suggests that there would exist an intrinsic relationship between the 
superconductivity and the As-bridged superexchange antiferromagnetic interactions. Here 
we have assumed that the value of J2 does not change with doping or pressure in comparison 
with the ones of the parent compounds. This can be verified from Fig. [TU Fig. [TH shows 
how the superexchange interactions J\ and J2 change with the pressure, from which we find 
that J 2 slightly changes with the pressure. However, J x drops quickly with increasing pres- 
sure. This is because the energy difference between the non-magnetic state and the square 
antiferromagnetic Neel state can be significantly reduced by the pressure, as shown in Fig. 

m 

In comparison with the bare and Jh-, Ji and J 2 are very small. When we study 
charge dynamics like charge transport happening in iron pnictides, t^-Jn part will thus 
play a dominant role. In Eq. ($Z§, if the Ji and J2 terms are ignore, the Hamiltonian 
will look similar to the double exchange model that was proposed for describing physical 
properties of colossal magneto-resistance (CMR) of manganese oxidesj^, l^]. However, 
there is an essential difference. In manganese oxides the crystal field splitting between e g 
and t 2g orbitals is very large, the moments polarized by the Hund's rule coupling on the t 2g 
orbitals are completely localized and only e g electrons can hop on the lattice. In contrast, in 
iron pnictides, all five 3d orbitals of Fe have contributions to the moment and at the same 
time they can also hop on the lattice. Nevertheless, we believe that their charge dynamics 
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FIG. 13: Various critical superconductivity temperature T c versus exchange interaction^ Jo, f° r 
AFe^Aso (A=Ba, Ca, or Sr) with doping or under high pressures. Tc are taken from Ref. 
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FIG. 14: Variations of superexchange antiferromagnetic interaction constants J\ and J2 with the 
imposed pressure for BaFe2As2- 

shares a common feature. In manganese oxides there is a strong spin-dependent scattering 
on the conduction electrons due to the strong Hund's coupling between e g and t 2g levels of 
Mn ions. If the localized core spins of Mn are aligned ferromagnetically, the scattering of the 
electrons due to the core spins will be dramatically reduced because of no spin-flip, leading 
to a giant magneto-resistance and a sharp drop of resistivity with decreasing temperature in 
the ferromagnetic phase. In iron pnictides, there is also a strong spin-dependent scattering 
on conduction electrons caused by the Hund's rule coupling. In low temperatures, the 
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moments of Fe ions are in the collinear antiferromagnetic order, in which the Fe moments 
are aligned antiferromagnetically along one direction but ferromagnetically along the other 
direction perpendicular. Similar as in manganese oxides, the scattering of electrons along 
the ferromagnetic direction is significantly reduced, leading to a sharp drop of resistivity 
and a large magneto-resistivity in the collinear antiferromagnetic phase, in agreement with 



experimental measurements 
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V. CONCLUSION 



In conclusion, we have reported calculated results on the electronic band structures of 
AFe2As2 (A=Ba, Ca, or Sr) by using the first-principles electronic structure calculations. 
The ground state of AFe2As2 is shown to be a collinear antiferromagnetic semimetal with a 
large magnetic moment around each Fe ion. The electronic structure is weakly affected by 
the spin-orbit interaction, but strongly altered by pressure. We have determined the density 
of states at the Fermi level, the specific heat coefficient and the Pauli susceptibility in both 
the non-magnetic and collinear antiferromagnetic states. The effective antiferromagnetic 
coupling constants of the Fe moments are also estimated assuming that the low-energy spin 
dynamics is approximately described by the Heisenberg model with the nearest and the 
next-nearest neighboring exchange terms. Based on the analysis of electronic and magnetic 
structures, we proposed that the low-energy physics of AFe 2 As 2 can be effectively described 
by the t — Jh — Ji — J2 model, defined by Eq. ([2]). 

This work is partially supported by National Natural Science Foundation of China and 
by National Program for Basic Research of MOST, China. 
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